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Abstract. The eclipse mapping method is an inversion technique that makes use of the information contained in 
eclipse light curves to probe the structure, the spectrum and the time evolution of accretion discs. In this review 
I present the basics of the method and discuss its different implementations. I summarize the most important 
results obtained to date and discuss how they have helped to improve our understanding of accretion physics, from 
testing the theoretical radial brightness temperature distribution and measuring mass accretion rates to showing 
the evolution of the structure of a dwarf novae disc through its outburst cycle, from isolating the spectrum of a 
D ■ disc wind to revealing the geometry of disc spiral shocks. I end with an outline of the future prospects. 

00 



OS 



> 
(N 



o 
o 



Or 



>< 



1 Introduction 



Accretion discs are cosmic devices that allow matter to efficiently accrete over a compact source by remov- 
ing its angular momentum via viscous stresses while transforming the liberated gravitational potential 
energy into heat and, thereafter, radiation fllTf . They are widespread in astrophysical environments, from 
sheltering the birth of stars to providing the energetics of quasars and active galactic nuclei. It is, however, 
CK ' in mass-exchanging binaries such as non-magnetic Cataclysmic Variables (CVs) that the best environ- 

ment for studies of accretion discs are possibly found. In these close binaries mass is fed to a white dwarf 
by a Roche lobe filling companion star (the secondary) via an accretion disc, which usually dominates 
the ultraviolet and optical light of the system [p7j . 

Accretion discs in CVs cover a wide range of accretion rates, M, and viscosity regimes. For example, the 
sub-class of dwarf novae comprises low- mass transfer CVs showing recurrent outbursts (of 2-5 magnitudes, 



on timescales of weeks to months) which reflect changes in the structure of the discs - from a cool, 
optically thin, low viscosity state to a hot, optically thick, high viscosity state - and which are usually 
parameterized as a large change in the mass accretion rate (M= 10" 11 Mq yr^ 1 i— ► 10~ 9 M© yr~ x ) |39f. 
On the other hand, nova-like variables seem to be permanently in a high viscosity state, presumably as 
a result of the fact that the accretion rate is always high. 

The temperatures in CV discs may vary from 5000 K in the outer regions to over 50000 K close to 
the disc centre, and the surface density may vary by equally significant amounts over the disc surface. 
Therefore, the spectrum emitted by different regions of the accretion disc may be very distinct. Addition- 
ally, the bright spot (formed by the impact of the gas stream from the inner Lagrangian point on the disc 
rim), the white dwarf at disc centre, and the secondary star may all contribute to the integrated light 
of the binary. Because what one directly observes is the combination of the spectra emitted from these 
diverse regions and sources, the interpretation of disc observations is usually plagued by the ambiguity 
associated with composite spectra. The most effective way to overcome these difficulties is with spatially 
resolved data. 

Two complementary indirect imaging techniques were developed in the 1980's that provide spatially 
resolved observational constraints on accretion discs on angular scales of micro arcseconds - well beyond 
the capabilities of current direct imaging techniques. One is Doppler Tomography p2|, which is treated in 
detail elsewhere is this book |51|. It uses the changes in line emission profile with orbital phase to probe 
the dynamics of accretion discs and is applicable to binaries over a large range of orbital inclinations, 
although it is restricted to emission line data. 

The other is Eclipse Mapping J21J . It assembles the information contained in the shape of the eclipse 
into a map of the accretion disc surface brightness distribution. While its application is restricted to 
deeply eclipsing binaries, eclipse mapping can be used with continuum as well as line data. When applied 
to time-resolved spectroscopy through eclipses this technique delivers the spectrum of the disc at any 
position on its surface. Information on the radial dependence of the temperature and vertical temperature 
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gradients (for optically thick regions), or temperature, surface density and optical depth (where the disc 
is optically thin) can be obtained by comparing such spectra with the predictions of models of the vertical 
disc structure. The spatial structure of the emission line regions over the disc can be similarly mapped 
from data of high spectral resolution. Moreover, an eclipse map yields a snapshot of an accretion disc at 
a given time. By eclipse mapping an accretion disc at different epochs it is possible to follow the secular 
changes of its radial brightness temperature distribution - for example, thought the outburst cycle of a 
dwarf nova - allowing crucial tests of accretion disc instability and viscosity models. 

After more than a decade of experiments, eclipse mapping has now become a mature and well estab- 
lished technique. There are already many good reviews on this topic in the literature |2^ , ^ , pTy6C[| . The 
main aim of this review is therefore not to provide another description of the technique, but to make a 
summary of the results obtained so far giving emphasis on the impact they had on our understanding of 
the physics of accretion discs. 

2 The maximum entropy eclipse mapping method 

2.1 The principles 

The three basic assumptions of the standard eclipse mapping method are: (i) the surface of the secondary 
star is given by its Roche equipotential, (ii) the brightness distribution is constrained to the orbital plane, 
and (iii) the emitted radiation is independent of the orbital phase. The first assumption seems reasonably 
robust. The others are simplifications that do not hold in all situations. A discussion on the departures 



from assumptions (ii) and (iii) is presented in section 2.5. 

A grid of intensities centred on the white dwarf, the eclipse map, is defined in the orbital plane. 
The eclipse geometry is specified by the inclination i, the binary mass ratio q (=Mi/M\, where M% and 
Mi are the masses of, respectively, the secondary star and the white dwarf) and the phase of inferior 
conjunction <fio f2UE3j. Given the geometry, a model eclipse light curve can be calculated for any assumed 
brightness distribution in the eclipse map. A computer code then iteratively adjusts the intensities in 
the map (treated as independent parameters) to find the brightness distribution the model light curve 
of which fits the data eclipse light curve within the uncertainties. The quality of the fit is checked with 
a consistency statistics, usually the reduced x 2 . Because the one-dimensional data light curve cannot 
fully constrain a two-dimensional map, additional freedom remains to optimize some map property. A 
maximum entropy (MEM) procedure |18|,[49| is used to select, among all possible solutions, the one that 
maximizes the entropy of the eclipse map with respect to a smooth default map. 

Figure [j] illustrates the simulation of the eclipse of a fitted brightness distribution while showing the 
comparison between the resulting model light curve and the data light curve. The geometry in this case 
is q = 0.3 and i = 81°. The left-hand panels show the data light curve (small dots) and the model 
light curve (solid line) as it is being drawn at five different orbital phases along the eclipse (indicated in 
the lower right corner). The right-hand panels depict the corresponding geometry of the binary for each 
orbital phase, in which the secondary star progressively occults the accretion disc as well as the white 
dwarf and the bright spot. The middle panels show the best-fit disc brightness distribution and how it is 
progressively covered by the dark shadow of the secondary star during the eclipse. At phase <f> = —0.08 
only a small fraction of the outer, faint disc regions are eclipsed and there is only a small reduction in 
flux in the light curve. The eclipse of the bright spot at the edge of the disc and of the bright inner disc 
regions occur at about the same time (slightly after <j> = —0.04) and coincide with the steepest ingress 
in the light curve. The flux at phase = does not go to zero because a significant fraction of the disc 
remains visible at mid-eclipse. The asymmetry in the egress shoulder of the light curve maps into an 
enhanced brightness emission in the trailing side of the disc (the right hemisphere of the eclipse map in 
Fig J). 

2.2 The expressions 

The expressions governing the eclipse mapping problem are as follows. One usually adopts the distance 
from the disc centre to the internal Lagrangian point, Rli, as the length scale. With this definition the 
primary lobe has about the same size and form for any reasonable value of the mass ratio q |21f| . If the 
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Fig. 1. Simulation of a disc eclipse (q — 0.3, i = 78°). Left-hand panels: data light curve (dots) and model light 
curve (solid line) for five different orbital phases (indicated in the lower right corner). Middle panels: eclipse maps 
in a false color blackbody logarithmic scale. Roche lobes for q = 0.3 are shown as dotted lines; crosses mark the 
centre of the disc. The secondary is below each panel and the stars rotate counter-clockwise. Right-hand panels: 
the corresponding geometry of the binary for each orbital phase. 



eclipse map is an N points flat, square grid of side XRli, each of its surface element (pixel) has an area 
(XRli) 2 /N and an associated intensity Ij. The solid angle comprised by each pixel as seen from the earth 
is then 
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where d is the distance to the system. The value of A defines the area of the eclipse map while the choice 
of N sets its spatial resolution. 

The entropy of the eclipse map p with respect to the default map q is defined as 
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where p and q are written as 

Pj = ^T ' <& = vwT' (3) 

Other functional forms for the entropy appear in the literature [£3|,|37|]. These are equivalent to (||) when 
p and q are written in terms of proportions. 

The default map Dj is generally defined as a weighted average of the intensities in the eclipse map, 



J2k ^Jk 



D j = ^ Jfc fc , (4) 



where the weight function ujjk is specified by the user. A priori information about the disc (p. ex., axi- 
symmetry) is i ncluded in the default map via ujjk- Prescriptions for the weight function ujjk are discussed 



in section 2^3. In the absence of any constraints on Ij, the entropy has a maximum S max = when 
Pj = qj, or when the eclipse map and the default map are identical. 

The model eclipse light curve ro((/>) is derived from the intensities in the eclipse map, 

N 

m(</>) = e 2 J2 IM4>) , (5) 

where <fi is the orbital phase. The occultation function Vj ((f)) specifies the fractional visibility of each pixel 
as a function of orbital phase and may include fore-shortening and limb darkening factors p3|I3|3(| . The 
fractional visibility of a given pixel may be obtained by dividing the pixel into smaller tiles and evaluating 
the Roche potential along the line of sight for each tile to see if the potential falls below the value of 
the equipotential that defines the Roche surface. If so, the tile is occulted. The fractional visibility of the 
pixel is then the sum of the visible tiles divided by the number of tiles. 

The consistency of an eclipse map may be checked using the \ 2 as a constraint function, 
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where d(4>) is the data light curve, <t(^>) are the corresponding uncertainties, r((jj) is the residual at the 
orbital phase (/>, and M is the number of data points in the light curve. Alternatively, the constraint 
function may be a combination of the x 2 and the R-statistics || , 
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to minimize the presence of correlated residuals in the model light curve J7] . For the case of uncorrelated 
normally distributed residuals, the R-statistics has a Gaussian probability distribution function with 
average zero and unity standard deviation. Requiring the code to achieve an R — 0, is equivalent to 
asking for a solution with uncorrelated residuals in the model light curve. 

The final MEM solution is the eclipse map that is as close as possible to its default map as allowed by 



the constraint imposed by the light curve and its associated uncertainties |23 43]. In matematical terms, 
the problem is one of constrained maximization, where the function to maximize is the entropy and 
the constraint is a consistency statistics that measures the quality of the fitted model to the data light 
curve. Different codes exist to solve this problem. Many of the eclipse mapping codes are based on the 
commercial optimization package MEMSYS Eq|. Alternative implementations using conjugate-gradients 
algorithms [0J§], CLEAN- like algorithms |5l] and, more recently, genetic algorithms |13| are also being 
used. 

2.3 Default maps 

A crucial aspects of eclipse mapping is the selection of the weight function for the default map, which 
allows the investigator to steers the MEM solution towards a determined type of disc map. A list of 
different prescriptions for uijk is given in Table [j]. 
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Table 1. Prescriptions for default maps 



weight functions 



reference 



A) most uniform map: uijk — 1 

B) smoothest map: uijk = exp 

C) most axisymmetric map: ujjk — exp 
(full azimuthal smearing) 

D) limited azimuthal smearing: uijk — exp 
(constant angle 8) 

E) limited azimuthal smearing: ujjk — exp 
(constant arc length s) 
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NOTES: djfc is the distance between pixels j and fc; Rj and -Rfc are the distances 
from pixels j and k to the centre of the disc; Ojk is the azimuthal angle between 
pixels j and fc; and Sjfc is the arc-length between pixels j and k. 



Choosing ujjk — 1 (option A) results in a uniform default map and will lead to the most uniform 
eclipse map consistent with the data. This happens not to be a good choice for eclipse mapping because 
it results in a map severely distorted by criss-crossed artifacts pd|pi]j5l| ] . This effect may be reduced by 
setting the weight function as a Gaussian profile of width A (option B), which results in the smoothest 
map that fits the data. 

The third case (option C) sets Dj as an axi-symmetric average of the eclipse map and will lead to the 
most nearly axi-symmetric map that fits the data. It suppresses the azimuthal information in the default 
map while keeping the radial structure of Ij on scales greater than Ar. This seems a reasonable choice for 
accretion disc mapping because one expects the disc material to be roughly in Keplerian orbits, so that 
local departures from axi-symmetry will tend to be diminished by the strong shear. This is a commonly 
used option and is also known as the default map of full azimuthal smearing. 

The full azimuthal smearing default results in rather distorted reproduction of asymmetric structures 
such as a bright spot at the disc rim. In this case, the reconstructed map exhibits a lower integrated flux 
in the asymmetric source region, the excess being redistributed as a concentric annulus about the same 
radial distance. By limiting the amount of azimuthal smearing it is possible to alleviate this effect and 
to start recovering azimuthal information in the accretion disc. Two prescriptions in this regard were 
proposed. Reference [J44| limited the amount of azimuthal smearing by averaging over a polar Gaussian 
weight function of constant angles along the map (option D) while |9l chosen to use a polar Gaussian 
function of constant arc length through the map (option E) . 

Figure H shows the effects of the three last weight functions when applied to an artificial map containing 
three Gaussian spots at different radial distances from its centre. The default with constant angles is more 
efficient to reproduce asymmetries in the inner disc regions (such as an accretion column, or the expected 
dipole pattern for velocity-resolved line emission mapping, or in presence of a bipolar wind emanating 
from the inner disc), while the default of constant arc length is more efficient in recovering asymmetries 
in the outer parts of the disc (such as a bright spot at disc rim, or in the case of an eccentric disc) . The 
choice between these two default functions, in a given case, is defined by whether it is more important 
to have a better azimuthal resolution at the inner or outer disc regions. The radial profile is not affected 
by this choice. 

Other possibilities concerning the default map were proposed in terms of the combination of different 
weight functions |ll|j5lj| . Particularly, the mix of the smoothest and the most axi-symmetric defaults is, 
in a sense, equivalent to the default of limited azimuthal smearing and leads to similar results. Another 
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Fig. 2. Effects of the different weight functions for the default map. (a) The original map, with three Gaussian 
spots. The corresponding default map obtained using the default of (b) full azimuthal smearing; (c) constant 
angles; and (d) constant arc length. From M. 

interesting proposal is that of a negative weight function, that may be used to avoid or minimize a certain 
map property (p. ex., the presence of the undesired criss-crossed arcs) |5jJ . 

2.4 The uneclipsed component 

Reference |17j found that the entropy function can be a useful tool to signal and to isolate the fraction of 
the total light which is not coming from the accretion disc plane. They noted that when the light curve 
is contaminated by the presence of additional light (e.g., from the secondary star) the reconstructed map 
shows a spurious structure in the regions farthest away from the secondary star (the upper lune of the 
eclipse maps in Fig. [l], hereafter called the 'back' side of the disc). This is because the eclipse mapping 
method assumes that all the light is coming from the accretion disc, in which case the eclipse depth and 
width are correlated in the sense that a steeper shape corresponds to a deeper eclipse. The addition of 
an uneclipsed component in the light curve (i.e., light from a source other than the accretion disc) ruins 
this correlation. To account for the extra amount of light at mid-eclipse and to preserve the brightness 
distribution derived from the eclipse shape the algorithm inserts the additional light in the region of the 
map which is least affected by the eclipse, leading to a spurious front-back disc brightness asymmetry. 
Since the entropy measures the amount of structure in the map, the presence of these spurious structures 
is nagged with lower entropy values. 

The correct offset level may be found by comparing a set of maps obtained with different offsets and 
selecting the one with highest entropy. Alternatively, the value of the zero-intensity level can be included 
in the mapping algorithm as an additional free parameter to be fitted along with the intensity map in 
the search for the MEM solution [ff|]45|| . A detailed discussion on the reliability and consistency of the 
estimation of the uneclipsed component can be found in [|9| . 

2.5 Beyond the standard assumptions 

The standard eclipse mapping assumes a simple flat, geometrically thin disc model. Real discs may 
however violate this assumption in the limit of high M. Disc opening angles of a ^ 4° are predicted 
for M^,5x 10 -9 Mq yr^ 1 |3^,|5(J. At large inclinations (i £ 80°) this may lead to artificial front-back 
asymmetries in the eclipse map (similar to those discussed in section E^J) because of the different effective 
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Fig. 3. Examples of reconstructed radial brightness temperature distributions with the flat-disc eclipse mapping 
for the case of flared discs. The inclinations i and disc opening angle a are indicated for each case. From [|43|. 



areas of surface elements in the front and back sides of a flared disc as seen by an observer on Earth. In 
extreme cases, this may lead to obscuration of the inner disc regions by the thick disc rim (e.g., |27j| ). 

Motivated by the front-back asymmetry that appeared in the flat-disc map and by the difficulties in 
removing the asymmetry with the assumption of an uneclipsed component, |40,4l| introduced a flared 
disc in their eclipse mapping of ultraviolet light curves of the dwarf nova Z Cha at outburst maximum. 
They found that the asymmetry vanishes and the disc is mostly axi-symmetric for a disc opening angle 
of a = 6°. 

Simulations JJ3[ show that eclipse mapping reconstructions obtained with the flat-disc assumption 
result in good reproduction of the radial temperature distribution of flared accretion discs provided that 
the inner disc regions are not obscured by the disc rim (Fig. 0). 

Reference pO| shows that it is usually impossible to distinguish between a flared disc and an uneclipsed 
component to the total light. Both effects lead to the appearance of spurious structures in the back regions 
of the disc, and eclipse maps obtained with either model may lead to equally good fits to the data light 
curve. Reference [0 pointed out that spectral eclipse mapping could help in evaluating the importance of 
each of these effects in a given case. If the uneclipsed component is caused by an optically-thin, vertically- 
extended disc wind, the uneclipsed spectrum shows a Balmer jump in emission plus strong emission lines, 
while in the case of a flared disc the spurious uneclipsed spectrum should reflect the difference between 
the disc spectrum of the back (deeper atmospheric layers seen at lower effective inclinations) and the 
front (upper atmospheric layers seen at grazing incidence) sides and should mainly consist of continuum 
emission filled with absorption lines. 

Because of the assumption that the emitted radiation is independent of the orbital phase, in the 
standard eclipse mapping method all variations in the eclipse light curve are interpreted as being caused 
by the changing occultation of the emitting region by the secondary star. Thus, out-of-eclipse brightness 
changes (p. ex., orbital modulation due to anisotropic emission from the bright spot) has to be removed 
before the light curves can be analyzed. The usual approach is to interpolate the out-of-eclipse light 
curve across the eclipse phases [ pip7[ . An alternative approach is to apply a light curve decomposition 
technique to separate the contributions of the white dwarf, bright spot and accretion disc B3,p3[ . This 
technique however requires high signal-to-noise light curves and good knowledge of the contact phases of 
the white dwarf and bright spot, which limits its application to a few objects. 

A step to overcome these limitations was done by |0] with the inclusion of a disc rim in the eclipse 
mapping method. The out-of-eclipse modulation is modeled as the fore-shortening of an azimuthally- 
dependent brightness distribution in the disc rim. This procedure allows to recover the azimuthal (phase) 
dependency of the bright spot emission. It however requires a good estimate of the outer disc radius. 
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Fig. 4. Reconstructing asymmetric brightness distributions with light curves of low signal-to-noise and incomplete 
phase coverage. The left-hand panel shows the artificial light curves (dots with error bars) and corresponding 
eclipse mapping models (solid lines). Horizontal dashed lines indicate the true zero level in each case. Vertical 
dotted lines mark ingress/egress phases of the white dwarf and mid-eclipse. The middle and right-hand panels 
show, respectively, the original maps and the reconstructions in a logarithmic greyscale. Bright regions are dark; 
faint regions are white. A cross marks the center of the disc; dotted lines show the Roche lobe and dashed lines 
depict the locus of the far edge of the shadow of the secondary star along the eclipse. The secondary is to the 
right of each map and the stars rotate counter-clockwise. 

The more advanced code of |42]| : including a flared disc, the disc rim, and the surface of the Roche-lobe 
filling secondary star, expanded the eclipse mapping method into a three-dimensional mapping technique. 
Nevertheless, it comes along with a significant increase in the degrees of freedom that aggravates the 
problem of non-uniqueness of solutions. 



2.6 Performance under extreme conditions 

In an eclipse mapping reconstruction, the brightness of a given surface element is derived from the 
information given by the changes in flux caused by its occultation (at ingress) and reappearance (at 
egress) from behind the secondary star. In the case of an eclipse light curve with incomplete phase 
coverage, there are regions in the disc for which only one of these pieces of information is available. 
Moreover, for a system with low inclination, there are regions in the back side of the disc which are never 
covered by the shadow of the secondary star and, therefore, there is no information about the brightness 
distribution of these regions on the shape of the light curve. This section presents simulations aiming to 
assess the reliability of eclipse mapping reconstructions obtained under the combined extreme conditions 
of incomplete eclipse coverage, low binary inclination and relatively low signal-to-noise data. 

Four artificial brightness distributions with asymmetric polar Gaussian spots on the back, front, 
leading and trailing sides of the disc were constructed (Fig. 0). A low- inclination geometry (q = 1 and 
i = 71°) was adopted to simulate the eclipses and light curves with signal-to-noise S/N ~ 5 — 15 and an 
incomplete set of orbital phases were produced. The artificial light curves were analyzed with the eclipse 
mapping method and the results are shown in Fig. 0. 

For the adopted set of orbital phases, the leading side of the disc (the lower hemisphere of the eclipse 
maps in Fig. |lj) is mapped by the moving shadow of the secondary star both during ingress and egress, 
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whereas much of the trailing side of the disc is only mapped by the secondary star at ingress phases. 
Dashed lines in the eclipse maps mark the locus of the far edge of the shadow of the secondary star along 
the eclipse. Regions to the left of this line are never covered by the secondary star. 

Despite the incomplete eclipse coverage, good quality reconstructions are obtained for the front, 
trailing and leading maps. The spots appear spread in radius due to the low signal-to-noise of the light 
curves, and are elongated in azimuth because of the intrinsic azimuthal smearing effect of the eclipse 
mapping method. The results are equally good for the leading and the trailing maps, despite the fact 
that the spot in the latter case is located in the disc region for which there is limited information in the 
shape of the light curve. For the back map, much of the asymmetric brightness distribution completely 
escapes eclipse. Not surprisingly, the eclipse map does not correctly reproduce the brightness distribution 
in the disc regions beyond those covered by the secondary star. The missing flux appears in the uneclipsed 
component. 

These simulations show that eclipse mapping obviously fails to recover the brightness distribution of 
disc regions for which there is no information in the shape of the eclipse, but performs reasonably well 
in the case of incomplete phase coverage even with relatively low signal-to-noise data. 

The brightness distribution of the back map approximately simulates the intrinsic front-back asym- 
metry of a flared disc as seen at a high inclination angle (i > 80°). It is fortunately that, in these cases, 
the shadow of the secondary star maps most (if not all) of the primary Roche lobe for any reasonable 
mass ratio. 

Simulations of reconstructions from light curves of more limited phase coverage are presented in H . 
Tests on the reliability of eclipse mapping reconstructions under a variety of other conditions can be 
found in the literature |,|,|ll|,[4|2l|,||,|l) . 

3 A summary of results 

When eclipse mapping appeared in the mid- 1 980 's, the standard picture of an accretion disc in a CV 
was that of a flat, nearly axi-symmetric disc with a bright spot on its edge. This section reviews some of 
the eclipse mapping results that helped to improve this picture either by allowing key tests of theoretical 
expectations or by revealing new and unexpected aspects of the physics of accretion discs. 

3.1 Classical results 

Early applications of the technique were useful to show that accretion discs in outbursting dwarf novae pij 
and in long-period novalike variables p9,E7J closely follow the expected radial dependence of temperature 
with radius for a steady-state disc, T ex i? -3 ' 4 , and to reveal that the radial temperature profile is 
essentially flat in the short period quiescent dwarf novae |6^|33|,Q (Fig. ||). This suggests that the 
viscosity in these short period systems is much lower in quiescence than in outburst, lending support to 
the disc instability model, and that their quiescent discs are far from being in a steady-state. Eclipse 
mapping studies also contributed to the puzzle about the SW Sex stars - a group of mostly eclipsing 
novalike variables with periods in the range 3-4 hr that display a number of unexplained phenomena - 
by showing that the radial temperature profile in these systems is noticeably flatter than the T ex R~ 3 / 4 
law ||,[47|]. A n& t radial temperature distribution is also suggested for the old novae V Per, that lies in 
the middle of the CV period gap |6l| . Unfortunately, this result is rather uncertain because it is based on 
uncalibrated white-light data and there is a large uncertainty on the eclipse geometry. A comprehensive 
account of these pioneering results can be found in J23| . 

It has been a usual practice to convert the intensities in the eclipse maps to blackbody brightness 
temperatures and then compare them to the radial run of the effective temperature predicted by steady 
state, optically thick disc models. A criticism about this procedure is that a monochromatic blackbody 
brightness temperature may not always be a proper estimate of the disc effective temperature. As pointed 
out by B, a relation between these two quantities is non-trivial, and can only be properly obtained by 
constructing self-consistent models of the vertical structure of the disc. Nevertheless, the brightness 
temperature should be close to the effective temperature for the optically thick disc regions. 

From the T(R) diagram it is possible to obtain an independent estimate of the disc mass accretion rate. 
Reference |23| compiled the inferred mass accretion rates from a dozen of eclipse mapping experiments 
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Fig. 5. The radial temperature profile of the dwarf nova Z Cha in outburst and in quiescence. The solid lines 
show steady-state blackbody disc models for mass accretion rates from 10 -8 to 10~ 12 M@ yr -1 . From p3|. 

to construct an MxP or j, diagram. An updated version of this diagram is shown in Fig. o. It seems a bit 
disappointing that the diagram is still loosely populated. In particular, there is yet no eclipse mapping 
estimate of M for a system inside the 2-3 hr CV period gap. There is a significant scatter in the M 
derived from different experiments for a given object (e.g., in UX UMa, from 1CP 8 ' 1 to 10 -8 ' A/qvi-" 1 
at 0.1 Rli)- Whereas part of this scatter is possibly a real effect due to long-term changes in the mass 
transfer rates, it stands as a warning that one should be careful in interpreting mass accretion rates 
derived from the brightness temperature distributions, as discussed above. 
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Fig. 6. Mass transfer rates at radii of 0.1 Rli (crosses) and 0.3 Rli (circles) as a function of the binary period. 



According to current evolutionary scenarios, CVs should evolve towards shorter orbital periods with 
decreasing mass transfer rates as a consequence of orbital angular momentum losses due mainly to 
magnetic braking (for systems above the period gap) or gravitational radiation (for systems below the 
gap) [ p7|j3q | . In Fig. o it appears that there is a tendency among the steady-state discs of novalike variables 
to show larger M for longer binary period - in agreement with the above expectation - and that the discs 
of novalike variables and outbursting dwarf novae have comparable M. The mass accretion rates in the 
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eclipse maps of novalike variables increase with disc radius. The departures from the steady-state disc 
model are more pronounced for the SW Sex stars (period range 3-4 hs) . Illumination of the outer disc 
regions by the inner disc or mass ejection in a wind from the inner disc are possible explanations for this 
effect. 

Multi-colour eclipse mapping is useful to probe the spectrum emitted by the different parts of the disc 
surface. Two-colour diagrams show that the inner disc regions of outbursting dwarf novae |lJ,|2J] and of 
novalike variables J|,|]J26| are optically thick with a vertical temperature gradient less steep than that of a 
stellar atmosphere, and that optically thin, chromospheric emission appears to be important in the outer 
disc regions (Fig. |7|). The fact that the emission from the inner disc regions is optically thick thermal 
radiation opens the possibility to use a colour-magnitude diagram to obtain independent estimates of the 
distance to the binary with a procedure similar to cluster main-sequence fitting. Distance estimates with 
this method were obtained for Z Cha §1, OY Car El, UU Aqr §, RW Tri || and UX UMa (§. 
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Fig. 7. Inferring the disc emission properties and the distance to Z Cha from the two-colour and colour-magnitude 
diagrams. Surface elements at the inner disc (R < 0.3 Rli) are represented by large dots, while elements in 
the outer disc regions are indicated by small dots. The solid and dashed curves in the right-hand panel show, 
respectively, the main sequence relationship for the best-fit distance and blackbody relationships for three different 
assumed distances. From [M. 



3.2 Spectral studies 

The eclipse mapping method advanced to the stage of delivering spatially-resolved spectra of accretion 
discs with its application to time-resolved eclipse spectrophotometry H] . The time-series of spectra is 
divided up into numerous spectral bins and light curves are extracted for each bin. The light curves are 
then analyzed to produce a series of monochromatic eclipse maps covering the whole spectrum. Finally, 
the maps are combined to obtain the spectrum for any region of interest on the disc surface. 

The spectral mapping analysis of the nova- like variables UX UMa J||3|,[44 45 1 and UU Aqr || shows 
that the inner accretion disc is characterized by a blue continuum filled with absorption bands and 
lines which cross over to emission with increasing disc radius (Fig. g|). The continuum emission becomes 
progressively fainter and redder as one moves outwards, reflecting the radial temperature gradient. Similar 
results were found for SW Sex p8| and RW Tri jlj|. Not surprisingly, these high-M discs seem hot and 
optically thick in their inner regions and cool and optically thin in their outer parts. 

However, the unprecedent combination of spatial and spectral resolution obtained with spectral map- 
ping started to reveal a multitude of unexpected details. In UU Aqr, the lines show clear P Cygni profiles 
at intermediate and large disc radii in an evidence of gas outflow M. In UX UMa, the comparison of 
spatially resolved spectra at different azimuths reveals a significant asymmetry in the disc emission at 
ultraviolet wavelengths, with the disc side closest to the secondary star showing pronounced absorption 
bands and a Balmer jump in absorption. This effect is reminiscent of that observed previously in OY Car, 
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Fig. 8. Spatially resolved spectra of the UXUMa accretion disc on August 1994 (gray) and November 1994 
(black). The spectra were computed for a set of concentric annular sections (mean radius indicated on the left, in 
units of Rli)- The most prominent line transitions are indicated by vertical dotted lines. From H. 



where the white dwarf emission seems veiled by an "iron curtain" p5j , and was attributed to absorption 
by cool circumstellar material || . The spectrum of the infalling gas stream in UX UMa and UU Aqr is 
noticeably different from the disc spectrum at the same radius suggesting the existence of gas stream 
"disk-skimming" overflow that can be seen down to R ~ (0.1 — 0.2) Rli- Spectra at the site of the bright 
spot suggest optically thick gas, with the Balmer jump and the Balmer lines in absorption. 

The spectrum of the uneclipsed component in these nova-like systems shows strong emission lines 
and the Balmer jump in emission indicating that the uneclipsed light has an important contribution 
from optically thin gas (Fig. 0). The lines and optically thin continuum emission are most probably 
emitted in a vertically extended disc chromosphere + wind |||(| ■ The uneclipsed spectrum of UX UMa 
at long wavelengths is dominated by a late-type spectrum that matches the expected contribution from 
the secondary star |J45| . Thus, the uneclipsed component seems to provide an unexpected but interesting 
way of assessing the spectrum of the secondary star in eclipsing CVs. 

Based on their spatially-resolved results, g suggested that the reason for the long standing discrep- 
ancies between the prediction of the standard disc model and observations of accretion discs in nova-like 
variables (e.g., pq , p4| , p6[ ) is not an inadequate treatment of radiative transfer in the disc atmosphere, 
but rather the presence of additional important sources of light in the system besides the accretion disc 
(e.g., optically thin continuum emission from a disc wind and possible absorption by circumstellar cool 
gas). 



3.3 Spatial studies 

Eclipse mapping has also been a valuable tool to reveal that real discs have more complex structures than 
in the simple axi-symmetric model. 

Besides the normal outbursts, short-period dwarf novae (SU UMa stars) exhibit superoutbursts in 
which superhumps develop with a period a few per cent longer than the binary orbital period. Normal 
superhumps appear early in the superoutburst and fade away by the end of the plateau phase. Late 
superhumps, displaced in phase by roughly 180° with respect to the normal superhumps, appear during 
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Fig. 9. The spectrum of the uneclipsed component in UX UMa (left) and UU Aqr (right, expressed as the 
fractional contribution to the total light). From |4aJa|. 



decline and persist into quiescence |57j). Eclipse mapping experiments have been fundamental in testing 
super hump models. 

Reference ||36f analyzed light curves of Z Cha during superoutburst with a modified eclipse mapping 
technique. Assuming that the superhump profile is fairly stable over a timescale of a dozen of binary orbits, 
he separated the eclipse of the superhump source by subtracting the light curve when the superhump 
maximum occurs far from eclipse from that in which the superhump is centred on the eclipse. The eclipse 
mapping of the resulting light curve show that the superhump light arises from the outer disc, and appears 
to be concentrated in the disc region closest to the secondary star. This result helped to establish the 
superhump model of |59[ , in which the normal superhumps are the result of an increased tidal heating 
effect caused by the alignment of the secondary star and an slowly precessing eccentric disc. 

Further evidence in favour of the existence of eccentric discs in CVs comes from the recent study of 
permanent superhumps in the short-period novalike variable V348 Pup by 142 ]. Their eclipse mapping 
analysis shows that the size of the disc emission region depends on superhump phase, and that the disc 
light centre is on the back side of the disc at superhump maximum. This phasing is reminiscent of that of 
the late superhumps in SU UMa stars. Their results indicate that the superhump maximum occurs when 
the secondary star is ligned up with the smallest disc radius, suggesting that these superhumps are the 
result of a modulation of the bright spot emission caused by the varying kinetic energy of the gas stream 
when it hits the disc edge J55fl . 

Ultraviolet observations of the dwarf nova OY Car in superoutburst show dips in the light curve 
coincident in phase with the optical superhump. Reference |1C| analyzed this data set with a modified 
version of the eclipse mapping method which simultaneously map the brightness distribution on the 
surface of the accretion disc and the vertical and azimuthal extent of a flaring at the edge of the disc. 
Their analysis indicates the presence of an opaque disc rim, the thickness of which depends on the disc 
azimuth and is large enough for the rim to obscure the centre of the disc at the dip phase (Fig. |10|) . These 
results are consistent with a model of normal superhumps as the consequence of time-dependent changes 
in the thickness of the edge of the disc, resulting in obscuration of the ultraviolet flux from the central 
regions and reprocessing of it into the optical part of the spectrum. Another evidence of discs with thick 
rims comes from the work of ]4l| , who found that a relatively large disc opening angle is required in order 
to explain the ultraviolet eclipse light curves of Z Cha in outburst. It seems that the discs of dwarf novae 
become flared during normal outbursts and that the thickening during superoutbursts may be sufficient 
for the disc rim to obscure the inner disc regions. 

Tidally induced spiral shocks are expected to appear in dwarf novae discs during outburst as the disc 
expands and its outer parts feel more effectively the gravitational attraction of the secondary star |3(J,|52| . 
Eclipse mapping of IP Peg during outburst H helped to constrain the location and to investigate the 
spatial structure of the spiral shocks found in Doppler tomograms p(J,p3| . The spiral shocks are seen in the 
continuum and C III+N III A4650 emission line maps as two asymmetric arcs of ~ 90 degrees in azimuth 
extending from intermediate to the outer disc regions (Fig. Kw. The Hell A4686 eclipse map also shows 
two asymmetric arcs diluted by a central brightness source. The central source probably corresponds to 
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Fig. 10. The dependency of the disc thickness with azimuth at three different wavelengths for OY Car in super- 
outburst. From 



the low-velocity component seen in the Doppler tomogram and is possibly related to gas outflow in a 
wind emanating from the inner parts of the disc. The comparison between the Doppler and eclipse maps 
reveal that the Keplerian velocities derived from the radial position of the shocks are systematically larger 
than those inferred from the Doppler tomography indicating that the gas in the spiral shocks has sub- 
Keplerian velocities. This experiment illustrates the power of combining the spatial information obtained 
from eclipse mapping with the information on the disc dynamics derived from Doppler tomography. 



3.4 Time-resolved studies 

Eclipse maps give snapshots of the accretion disc at a given time. Time-resolved eclipse mapping may be 
used to track changes in the disc structure, p. ex., to assess variations in mass accretion rate or to follow 
the evolution of the surface brightness distribution through a dwarf nova outburst cycle. 

The observed changes in the radial temperature distribution (and mass accretion rate) of eclipse maps 
obtained at different epochs in the high viscosity, steady-state discs of the novalike variables UX UMa 
and UU Aqr are evidence that the mass transfer rate in these system is variable 0|J] . 

Eclipse maps of Z Cha during superoutburst show a bright rim in the outer disc regions which decreases 
in brightness relative to the inner regions as the superoutburst proceeds and the superhumps fade awa y 
|p8| - This underscores the indications that the superhumps are sited at the outer disc rim (section |3^). 

Reference J46J obtained eclipse maps of the dwarf nova OY Car along the rise to a normal outburst. 
Their maps show that the outburst starts in the outer disc regions with the development of a bright ring, 
while the inner disc regions remain at constant brightness during the rise. The flat radial temperature 
profile of quiescence and early rise changes, within one day, into a steep distribution that matches a 
steady-state disc model for M= 10~ 9 M yr -1 at outburst maximum. Their results suggest that an 
uneclipsed component develops during the rise and contributes up to ~ 15 per cent of the total light at 
outburst maximum. This may indicate the development of a vertically-extended (and largely uneclipsed) 
disc wind, or that the disc is flared during outburst (see section 3.3). 

Time-resolved eclipse mapping covering the decline of an outburst and of a superoutburst were ob- 
tained, respectively, for IP Peg Q and OY Car 0|. In both cases the radial temperature distribution 
evolves with the inward traveling of a transition front that leaves behind a cool disc (T& ~ 5000 — 6000 K) 
while the temperatures at the inner disc remain almost constant at a higher value. The derived speed of 
this cooling front is ~ 0.14kms _1 for OY Car and ~ 0.8 kins" 1 for IP Peg. 

Eclipse maps covering the full outburst cycle of the long-period dwarf nova EX Dra JE| show the 
formation of a one-armed spiral structure in the disc at the early stages of the outburst |1| and reveal 
how the disc expands during the rise until it fills most of the primary Roche lobe at maximum light 
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Fig. 11. Eclipse mapping of spiral shocks in IP Peg. The light curves are shown in the left-hand panels and the 
eclipse maps are displayed in the middle and right-hand panels in a logarithmic greyscale. The notation is the 
same as in Fig. m From [H]. 

(Fig. [l2]). During the decline phase, the disc becomes progressively fainter until only a small bright region 
around the white dwarf is left at minimum light. The evolution of the radial brightness distribution 
suggests the presence of an inward and an outward-moving heating front during the rise and an inward- 
moving cooling front in the decline (Fig. |l3| ) . The inferred speed of the outward-moving heating front is 
of the order of 1 kms" 1 , while the speed of the cooling front is a fraction of that - in agreement with the 
results from IP Peg and OY Car. Their results also suggest a systematic deceleration of both the heating 
and the cooling fronts as they travel across the disc, in agreement with predictions of the disc instability 
model |33|]. A similar effect was seen in OY Car |TJ|. The radial temperature distributions shows that, as 
a general trend, the mass accretion rate in the outer regions is larger than in the inner disc on the rising 
branch, while the opposite holds during the decline branch. Most of the disc appears to be in steady-state 
at outburst maximum and, interestingly, also during quiescence. It may be that the mass transfer rate in 
EX Dra is sufficiently high to keep the inner disc regions in a permanent high viscosity, steady-state. An 
uncclipsed source of light was found in all maps, with a steady component associated to the secondary 
star and a variable component that is proportional to the out of eclipse brightness. Although disc flaring 
is likely in EX Dra during outburst, it seems it is not enough to account for the amplitude of the variation 
of the uneclipsed source. The variable component was therefore interpreted as emission arising from a 
disc wind, the strength of which depends on the disc mass accretion rate. 



4 Future prospects 

Eclipse mapping is a powerful probe of the radial and the vertical disc structures, as well as of the physical 
conditions in accretion discs. Partly due to the many experiments performed over the last 15 years, we 
have enriched our picture of accretion discs with an impressive set of new details such as gas outflow in 
disc winds, gas stream overflow, flared discs with azimuthal structure at their edge, ellipsoidal processing 
discs, sub-Keplerian spiral shocks, and moving transition fronts during disc outbursts. 
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Fig. 12. Sequence of eclipse maps of the dwarf nova EX Dra. The eclipse maps capture 'snapshots' of the disc 
brightness distribution in quiescence (h), on the rise to maximum (a-b), during maximum light (c), through the 
decline phase (d-f), and at the end of the eruption, when the system goes through a low brightness state before 
recovering its quiescent brightness level. The notation is the same as in Fig. m From pi. 



This is however far from being the end of the road. There are still many key eclipse mapping experi- 
ments remaining to be done. The spectral mapping of dwarf novae in outburst offers a unique opportunity 
to probe the physical conditions of disc spiral shocks and to critically test the disc instability model by 
comparing the spectra of disc regions ahead and behind the transition fronts. Eclipse mapping estimates 
of M for CVs inside the period gap may be instrumental in testing the current theories about the origin 
of the period gap. A still untouched area is the mapping of the flickering sources. Our understanding of 
this fundamental signature of accretion processes can certainly be considerably improved with flickering 
mapping experiments. The unprecedent combination of high spatial and high spectral resolution on the 
disc surface that can be achieved from time-resolved spectroscopic data yields an unbeatable amount of 
information to test and improve the current disc models. Fitting disc atmosphere models to the spatially- 
resolved spectra is the obvious next step to the spectral mapping experiments and will certainly be 
rewarding. 

Eclipse mapping have been continuously expanding into new domains. Of promising prospects are 
the direct mapping of physical parameters in the accretion disc [ p4[ , the mapping of the brightness 
distribution along the accretion stream and accretion column in magnetic CVs |29] , the 3-D mapping of 
flared discs and the surface of the secondary star p5p3] , and the combination of eclipse mapping and 
Doppler tomography 113] . 

The leap in knowledge about accretion physics that have been (and will probably continue to be) 
obtained by the study of close accreting binaries with tomographic techniques can lead to a better 
understanding of many other astronomical scenarios in which accretion discs may play an important role, 
such as AGNs and quasars - the environment of which are apparently more complicated and comparatively 
less well known than that of these compact binaries. 
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Fig. 13. Left: The radial intensity distributions of EX Dra through the outburst. Labels are the same as in Fig. [12. 
Dashed lines show the 1-a limit on the average intensity for a given radius. A dotted vertical line indicates the 
radial position of the bright spot in quiescence. Large vertical ticks mark the position of the outer edge of the 
disc and short vertical ticks indicate the radial position of a reference intensity level. Right: The radial brightness 
temperature distributions. Steady-state disc models for mass accretion rates of log M= —7.5, —8.0, —8.5, and 
—9.0 Mq yr _1 are plotted as dotted lines for comparison. The dot-dashed line marks the critical temperature 
above which the gas should remain in a steady, high mass accretion regime. The numbers in parenthesis indicate 
the time (in days) from the onset of the outburst. From M. 
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